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ABSTRACT 

The sensitivity of interferometers with linear polarizers to the CMB E and B mode are 
variant under the rotation of the polarizer frame, while interferometer with circular 
polarizers are equally sensitive to E and B mode. We present analytically and numeri- 
cally that the diagonal elements of window functions for CMB E/B power spectra are 
maximized in interferometric measurement of linear polarization, when the polarizer 
frame is in certain rotation from the associated baseline. We also present the simulated 
observation to show that the la errors on E/B mode power spectrum estimation are 
variant under the polarizer frame rotation in the case of linear polarizers, while they 
are invariant in the case of circular polarizers. Simulation of the configuration similar 
to the DASI shows that minimum la error on B mode in interferometer measure- 
ment with linear polarizers is 26 per cent of that in interferometric measurement with 
circular polarizers. The simulation also shows that the E/B mixing in interferometer 
measurement with linear polarizers can be as low as 23 per cent of that in interfer- 
ometric measurement with circular polarizers. It is not always possible to physically 
align the polarizer frame with all the associated baselines in the case of an interfer- 
ometer array (N>2). There exist certain linear combinations of visibilities, which are 
equivalent to visibilities of the optimal polarizer frame rotation. We present the linear 
combinations, which enables B mode optimization for an interferometer array (N>2). 

Key words: - cosmology: cosmic microwave background - techniques: interferometric 



1 INTRODUCTION 



The Cosmic Microwave Background (CMB) is expected 
to be linearly polarized by Thomson scattering at the 
last scattering surface and after re-ionization. The detec- 
tion of the CMB polarization has been rep orted by the 
Degree Angular Scale Interferometer (DASI) (|Kovac et al.l 
|2002| ) and recently by the Wilkinson Microw ave Anisotropy 
Probe (WMAP) satellite (|Page et al.ll2006t ). Measurements 
of the CMB temperature anisotropy with interferometers 
are made in the Very Small Array (VSA), the Cosmic Back- 
ground Imager (CBI) and many other experiments. The 
CMB polarization measurements with interferometers are 
on-going and planned in the experiments such as DASI, CBI 
and the Millimeter-wave Bolometric Int erferometer (MBI) 
jTucker et aHl2003l ; lKo"rotkov et al.1120061 ). With many desir- 
able features of an interferometer, interferometers are more 
and more employed in CMB polarization experiments. 

The CMB polarization can be decomposed 
into gradient-like E mode and curl-like B mode 
|Zaldarriaga and Seliakl Il997l ). B mode polarization is 
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not induced by scalar density perturbation but only tensor 
perturbation, while E mode p olarization is induced by both 
(|Seliak and Zaldarriagal Il997). Since tensor-to-scalar ratio 
is much smaller than unity in most inflationary models, B 
mode polarization is expected to be much smaller than E 
mode polarizatio n. Though there is com plication by grav- 
itational lensing (|Okamoto an d Hu 2003j), measurement of 
B mode polarization makes it possible pro bing the Univers e 
on the energy scale at inflationary period (Dodelson 2003). 

The la error on the para meter estimation can be fore- 
cast from the Fisher matrix (iDodelsonl 2003). It will be 
shown that in interferometric measurement of Stokes pa- 
rameter Q or U, the la error on E and B power spectra 
estimation varies with rotation of its polarizer frame. We 
will show that certain rotation of the polarizer frame from 
the associated baseline minimize the la error on either E 
or B mode power spectra estimation. For a feedhorn array 
(N>2), it is not always possible to realize the specific rota- 
tion of the polarizer frame from all the associated baselines. 
We will show that forming certain linear combinations of 
polarimetric visibilities is identical with physically rotating 
the polarizer frame to the optimal orientation, thereby en- 
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abling the B mode measurement optimization for a feedhorn 
array (N>2). 

This paper is organized as follows. We discuss Stokes 
parameters in §2. The formalism for interferometric CMB 
polarization measurement on a spherical sky is presented 
§3. In §4, with flat sky approximation, we show that inter- 
ferometric measurement of Stokes parameter Q or U can be 
configured to suppressing E mode, leading to smaller leakage 
between E and B mode. In §5, we discuss the variation of 
la error on power spectra estimation and show that B mode 
sensitivity is maximized, when the polarizer frame is in cer- 
tain rotation from the baseline. In §6, we show that there 
exists certain linear combinations of visibilities, which en- 
ables the B mode optimization for a feedhorn array (N>2). 
In §7, the summary and conclusion are given. In appendix, 
we show analytically that the diagonal element of window 
functions are maximized, when the polarizer frame is rotated 
from the baseline by certain angles. 



2 STOKES PARAMETERS 

There are St okes parame ters, which describe the state of 
polarizatio n (IKrauslll986|). which are me asured in reference 
to (eg,e$) (|Zaldarriaga and S cliak 1997). eg and are the 
unit vectors of the spheri cal coordinate system and given by 
jArfken and Weberlbood ) 

eg = i cos#cos</)+j cos 9 sin cj> — k sin#, 
= — i sin^+j cos0. 

Since Thomson scattering does not generate circular po- 
larization but linear polarization in early universe, the phase 
difference between Eg and E$ is zero and circular polariza- 
tion state V is redundant in the study on the CMB polar- 
ization. Stokes parameters Q and U are as follows: 

Ti2 t-,'2 \ 



Q = (E'l-El) 
U = (2E e E^) , 



(1) 

(2) 



where (. . .) indicates time average. Q and U transform un- 
der rotation of an angle ip on the plane perpendicular to 
direction n as 



Q'(n) = Q(n) cos2ip + U(n) sin 2^>, 
U'(n) = -Q(n)sin2i/> + [/(n)cos2^, 



(3) 
(4) 



with which the following q uantities can be constructed 
ijZaldarriaga and Scliak 1993): 



Q'(n) ± iU'(n) = e T ^(Q(n) ± iU(h)) 



(5) 



For all-sky analysis, Q an d U are expanded in terms o f 
spin ±2 spherical harmonics (Zaldarri aga and Seliak|[l997I ) 
as follows: 



Q{n)+iU{n) = 
Q(n)-if/(n) = 



I ,m 

£ 



-{aEjm + i a-B.im) iYi m (&), (6) 



) - 2 Y lm (h), (7) 



3 INTERFEROMETRIC MEASUREMENT 

The discussion in this section is for an ideal interferome- 
ter. An interferometer measures time-averaged correlation 



of two electric field from a pair of identical apertures posi- 
tioned at ri and at r2. The separation, B = ri — r-2, of two 
apertures is called the 'b aseline' and th e measured correla- 
tion is called 'visibility' (|Lawsonl Il999*) . Depending on the 
instrumental configuration, visibilities are associated with 
{El - Ey), (2E x E y ) and (E 2 X - E 2 y ±%2E X E V ) respectively, 
where x and y are axes of the polarizer frame. As discussed 
in fj2j Stokes parameters at angular coordinate (9,(f>) are de- 
fined in respect to two basis vectors eg and . Consider the 
polarization observation, whose antenna pointing is in the 
direction of angular coordinate (9a,<Pa)- The polarizers and 
baselines are assumed to be on the aperture plane. Then, the 
global frame coincides with the polarizer frame after Euler 
rotation (4>a, 9 a, VO on the global frame, where ip is the ro- 
tation around the axis in the direction of antenna pointing. 
Most of interferometer experiments for the CMB observation 
employ feedhorns for beam collection. After passing through 
the feedhorn system, an incoming off-axis ray becomes on- 
axis ray. Then, the basis vectors eg and e,/, of the ray after 
the feedhorn system are related to the basis vectors e x and 
e y of the polarizer frame as follows: 

" lV, (e flA + je^ A ) = e^* - ^^ + ie^), 



+ i e y = e 



where $ is given by 



$ = tan 



(8) 



sin 9 sin(0 — 4>a) 



sin a cos a a cos( 



6a) — cos 9 sin 9a 



+ tan 



sin 9a sin((p — <f>A) 



- sin 9 cos 9a + cos 9 sin ( 



Refer to Appendix [A] for the details on the derivation of 
With Eq. [8] we can easily show that 

(El - E 2 y ) + i(2E x E y ) = e-^~ 2 *\(E 2 g - E%) + t(2Eg E*)). 

With the employment of linear polarizers, the visibilities 
associated with (El — E 2 ) or (2E x E y ) are as follows: 

V Q , = f{v) J dOA(ft-fU) (9) 

xRe[e- l(2,/ '- 2 * (fl)) (Q(n) + i[/(n))] e 1 27ru ' fl , 
Vw = f(y) I dQA(n-n A ) (10) 



\Q(ii)+iU(n))], 



xlm [e 



where ua is the unit vector in the direction of antenna point- 
ing and f(v) is the frequency spectrum of the CMB polar- 
ization. □ With the employment of circular polarizers, the 
visibilities associated with (El — E^ ± i2E x E y ) are as fol- 
lows: 



(11) 



V RL = f(v) J d04(n-nA), 

x [Q(n) + l [/(n)]e l(2 ' ru,fl - 2 ^ +2 * (fl)) , 
Vlr = f(v) J d04(n-riA), 

X [Q(fi) - j f/ ( ft )] e *(2 7 ru.n+2V.-2<I.(n)) ) 

where R and L stand for right/left circular polarizers. 



(12) 



8B(v,T) I 

ST \t=T 
To is the CMB monopole temperature 



1 m 



where B(u, T) is the Plank function and 
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With Eq. [6] and [7] visibilities are expressed in terms of 
E/B mode as follows: 

V Q , = -f{y) J dQA(h-n A )e i( - 2nu A) (13) 

wT> r -s(2j*-2*(n)) I . ■ \ \f ~\ 

xRe|e yv (a E ,im + laajm) aximj , 

Vu' = -/M / dfM(n - r^) e l(2 ™' ft) (14) 

xlm [ e _ *( 2 ^ -2 *w) (a Bi!m + ia B ,i m ) 2 Vim] , 
Fbl = -/(f) fdfiA(n-nA) (15) 

w C i • \ \s i(27ru-ft-2;/; + 2$(ft)) 

X (CLE.lm + I ClB.l m ) 2 I lm& j 

Vis = -f{u) J AVLA{h-h A ) (16) 

/ % i(27ru-fl+2i/>-2$(fl)) 

x(a B , im - ia B ,im) -iYi m e v K ". 



4 REDUCING LEAKAGE BETWEEN E AND B 
MODE 

For the observation of small patch of sky, flat sky approxi- 
mation in small angle limit can be used. In flat sky approx- 
imation, Eg. l9l [101 1111 and 1121 are as follows: 



v Q , = m j dx 2 ^(x) 

xRe [e~ M *{Q(x) + i C/(x)}] e 12 ™'* , 
V v > = f(v) J dx 2 ,4(x) 

xlm [e- 2l *{Q(x) +W(x)}] e l2 ™' x , 
V RL = /M /dx 2 A(x)[Q(x)+iC/(x)]e l ^ u ' x - 2 ^, 



VLh = /(^) / dx 2 A(x)[Q(x)-i(7(x)] e l(2 " u ' x+2 ' / ' ) , 



where A(x) is the function of beam power pattern. Since 
visibilities are the convolution of the Fourier transform 
of a beam function with the Four ier transform of Q and 
[/ |M.P.Hobson and M aisinger 2002), they can be written as 
follows: 



V Q , = f[v) J d 2 u'A(u-u') 

x [cos(2 ip)Q(u) + sin(2 ip) U(u% 

V u , = f(u) /dVi(u-u') 

x [- sin(2 V)Q(u') + cos(2 ?/>) U(u% 
Vrl = f(u) J d 2 u'A(u-u')[e- l2 ^(Q(u') + iU(u'))], 

Vlh = f{y) /dVA(u-u') [e i2 *(Q(u')-il7(u'))], 



where tilde " indicates Fourier transform. In the flat sky ap- 
proximation in small angle li mit, Stokes parameter Q and U 
can be decomposed as follows (Zal darriaga and Seliakll9'97h : 

Q(u) = cos(20 u )S(u)-sin(2^u)B(u)), (17) 
!7(u) = sin(2^ u ) £(u) + cos(2(/> u ) B(u)), (18) 



where U is the direction angle of a vector u. With Eq. 1171 
and 1 181 visibilities are expressed in terms of E and B mode 
as follows: 

V Q , = f(u) J dVi(u-u') (19) 
x[cos(2(</> - M)E(u') + sin(2(V> - <p u ,)) B(u% 

V v , = f{y) J dVI(u - u') (20) 
x [- sin(2(</> - <t>u>))E(u) + cos(2(V> - U O) B(u% 

V RL = f[u) J d 2 u'i(u - u') (21) 
x [ e -* 2 ^-^'\E(u') + iB(u'))] , 

V LR = f(u) J d 2 u'i(u - u') (22) 

x [e ,2( *-*»' ) (£;(u')-«B(u'))] • 

A Gaussian beam is a good approximation to many CMB 
experiments and the Fourier transform of a Gaussian beam 
is exp[- |u ~ u 2 ' |2 ' T2 ], where a = 0.4245 FWHM. With such a 
beam, biggest contribution comes from u' = u in the inte- 
gration over u'. We can approximate visibilities as follows: 

Vq, « f{y) (23) 

x [cos(2(V> - 4>u))E{vl) + sin(2(V> - <j>*)) B(u)], 
V W « /(f) (24) 
x [- sin(2(V> - <pn))E(u) + cos(2(V> - <p u )) B(u)], 
Vrl « /(i/)[e- i2 <*-*^(£(u) + iB(u))], (25) 
Vb H « /(i/)[e° ( *- w (£(u)-iB(u))]. (26) 

As seen in Eq. [25] and 1261 the measurement of Vrl and 
Vlr measures E and B mode equally, independent of the 
polarizer rotation. From Eq. [23] and 1241 it is seen that Vq' 
and Vjji gets unequal contribution from E and B mode, when 
the baseline and the polarizer frame are aligned as follows: 



Vq, 

V v , 



f{v)E{n) : V = 0u 

f{v)B{u) : ^ = u + tt/4 

/(f) B(u) : V = 0u 

-/(i/)l(u) : </> = U + tt/4 



Physically, ip = U corresponds to aligning x axis of the 
polarizer frame with the baseline, and ip = <j> u + n/4 is ro- 
tating the a; axis of the polarizer frame from the baseline 
by 45° on the aperture plane. Compared with those of Vrl 
and Vrl in those configuration, either E or B mode in Vq, 
and Vu 1 is suppressed while the other mode is intact. The 
complete separation of E mode from B mode is not possible 
unless a full-sky m ap is made with infinite angular resolution 
l|Bunn et al.ll2003h . The leak from E mode into much weaker 
B mode causes serious problem. We can reduce E mode leak 
into B mode by measuring Vq, with ip = <f> u + tt/4 and Vu 1 
with ip = cj> u , in which E mode contribution is su ppressed. 

W e have computed the 2x2 leakage matrix l|Tegmarkl 
2001) for the simulated observation in !j5] Leb indicates 
the B mode leakage into E mode measurement and Lbe 
indicates the E mode leakage into B mode measurement. If 
Leb = Lbe = 0, there is no leakage at all l]Tegmarkll200ll ). 
Leb and Lbe are shown for various ip — U in Fig. [T] where 
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Figure 1. Lbb an d Lbe f° r various ip - 



Figure 2. Window functions in Vqi measurement with 3.4° 
FWHM 



Leb and Lbe are dimensionless, and visually identical. It 
shows that the leakage in Vqi (Vjji) is smallest at ip = </> u + 
7t/4 (ip = cp u ). It is seen that with the choice of ip — <p u , the 
leakage of Vqi and Vjji can be as low as 23 per cent of that 
of Vrl and Vlr- 



5 VARIATION OF la ERROR 

It has becom e standard to estimate band power spectra 
l|Knoxlll999f l by maximum likelihood method from CMB 
interferometric observations (|Hobson and Maisingerl [2002). 
With the Gaussianity of polarimetric visibility and noise, 
likelihood function is given by 



(2tt)-|S + N| 



exp[--A(S + N)- 1 A], 



where A is data, S is signal covariance matrix and N is noise 
covariance matrix. By exploring parameter space to maxi- 
mize likelihood, we can estimate parameters within certain 
error limits. The diagonal element of dS/dCi, where C; is 
the angular power spectrum, shows the sensitivity of the 
experiment over multipoles. Through this paper, a win- 
dow function means Wi = dS/dCi. In Appendix|B] we have 
shown analytically that the diagonal W t BB of Vqi is maxi- 
mized at ip — U ± 7r/4 while that of Vv> is maximized at 
Tp — 4> u , </> u ±7r/2. We have numerically computed Wj and 
W t BB for an interferometer of a 25 cm baseline with 30 ~ 
31 GHz signal frequency range and 3.4° Full Width at Half 
Maximum (FWHM) beam. They are shown in Fig. [2] and 
[3] and agree with the analytical result in Appendix [B] The 
window functions are normalized so that the peak value of 
the highest window function is a unit value. As also shown in 
Fig. [2] and [3] an interferometer is sensitive to the multipole 
range / ~ 2nu± Al/2, where u is a baselinelength divided by 
wavelength and Al is FWHM of the window function (for 
a circular Gaussian beam, A/ = 4V21n2/0 FW HM). In maxi- 
mum li kelihood est imation, it is usual to estimate the band 
powers (|Knoxlll999h . which are assumed to be flat over some 
multipole range. In an interferometer experiment, E and B 
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Figure 3. Window functions in Vjji measurement with 3.4° 
FWHM 



mode band power, Xee and Xbb are assumed to be flat over 
the multipole range the interferometer is sensitive to. 

The minimum possible variance on the par ameter es- 
timat ion can be forecast from the Fisher matrix (|Dodelsonl 
2003), which is defined as 



d 2 (ln£) 
dXidXj 



(27) 



N) 



(S + ISO- 



Evaluated at the maximum of the likelihood, the square 
root of diagonal element of the inverse Fisher matrix yields 
the marginalized la error on the parameter estimation. As 
shown analytically in appendix [B] and numerically in Fig. 
[2] and O diagonal windows functions of Vqi and Vw are 
maximized at certain ip — (p u . Therefore, AXee and AXbb 
are expected to be smallest with certain tp — (f> u . We have 
numerically computed AXee and A Ass from simulated ex- 
periments for various ip—(p u , which are shown in Fig.[4]and[5] 
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and at ip = cp u in Vy measurement. AAbb and A Abb from 
Vrl or Vlb. are invariant under the rotation of the polarizer 
frame. For the same assumed noise variance, the minimum 
A Abb from Vqi or Vjji is 26 per cent of AAbb from Vrl 
or Vlr- It is shown in Fig. [4] that AAbb is minimum at 
ip = <^u in Vqi measurement and at ip = U + 7r/4 in Vj// 
measurement. The minimum AAbb from Vqi or Vy/ is 1.2 
times bigger than the AAbb from Vrl or Vlb,. Vrl and Vrl 
have information on both of Q' and U', while Vqi and Vj// 
have information on either of Q' and U' . It may seem odd 
that AAbb of Vqi and Vj// can be smaller than AAbb of 
Vrl and Vrl, though more informations are contained in 
Vrl and Vrl- AAbb is the estimation error marginalized 
over Abe , since the likelihood in the simulated observation 
is the function of two parameters, Abe and Abb- AAbb is 
given by \° BB - X BB , where 



dC B {\i 



= 0, 



£b(-^bb) 



-1/2 



jCb(Xbb), which is the likelihood function marginalized over 
E mode band power, is as follows: 

X3b(Abb) = J CLXee C(\ee, Abb, Abb = 0). 

jC(Xee, Abb) becomes less sensitive to the variation of Abb 
with reduced contribution of E mode to visibilities. It makes 
the C B (\ BB ) more sharply peaked around \° BB , which leads 
to the reduction of AAbb. 

The lcr error on power spectra esti mation is the sum 
of sample variance and noise variance l|Park et al.l 120031 ) : 
AC; ~ Ci/(2l + 1) + N, where N is noise variance. In the 
case of circular polarizers, Abb is smaller than Abe due to 
the smaller sample variance of B mode than that of E mode, 
even though visibilities with circular polarizers have equal 
sensitivity to both E and B mode. 
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Figure 5. AAbb for various ip — <j> u 

In the simulation, we used the C ode for Anisotropi es in the 
Microwave Background (CAMBl lLewis et al. 2000) to com- 
pute the power spectra of ACDM with the tensor-to-scalar 
ratio (r = 0.3). The baseline length 25 [cm] is assumed with 
the signal frequency range 30~31GHz with 1GHz bandwidth 
so that the interferometer p robes the multipo le range of 
roughly El/Bl band of DASI (<Kovac et al.ll2002Tl. The noise 
covar iance matrix is assumed to be diagonal fWhite et al.l 
1 19991 ) and have a uniform value, for whi ch we assumed th e 
sensitivity of DASI: 60 x .1 2 tt Jy s 1/2 m 2 l|Prvke et al.ll2002h . 
We assumed the probe of three fields, three hundred sixty 
five days integration time for each field and simultaneous 
thirty six orientations for each baseline length. The equa- 
torial coordinates of the assumed three fields are (80°, 0°), 
(80°, 120°) and (80°, 240°). Fig. Hand [5] shows AAbb and 
AAbb for various ip — U . To be compared with the CAMB 
power spectra, AAbb and AAbb should be multiplied with 
1(1 + 1)/ (27r), where I = 162 is the multipole our assumed in- 
terferometer is most sensitive over. As shown in Fig. [5] that 
AAbb is minimum at ip = (p u + 7r/4 in Vqi measurement 



6 FEEDHORN ARRAY 

For optimization, the polarizer frame should be oriented 
with certain rotation from all the baselines formed by the 
feedhorn the polarizer is associated with. When there are 
just two feedhorns, there is a single baseline and orienting 
the polarizer frame to specific rotation physically is trivial. 
But in most of real CMB experiments, an array consisting of 
more than two feedhorns are employed. When all the feed- 
horns (apertures) are aligned on a straight line, orienting 
polarizers with specific rotation from all the associated base- 
lines is also trivial even for a (N>2) feedhorn array. When 
the feedhorn array is configured on so me two dimensional 
pattern for the optima l uv coverage (|Guvon and Roddierl 
l200ll ; |p~rvke et al"1l2002l ). it is not possible to realize specific 
rotation of the polarizer frame from all the associated base- 
lines. We can achieve the optimal polarizer rotation for a 
array of multiple feedhorns (N>2) by forming certain linear 
combination of visibilities as follows. With the employment 
of orthomode transducer (OMT), each baseline measures 
the real and imaginary part of Vqi and Vjy< respectively. 
Four measured values are obtained such that Vi = Re[VQ/], 
V 2 = lm[V Q i], V 3 = Re[Vu>], V 4 = lm\V V i}. With Eq.|l]and 
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I lOt Vi, V2, Vis, and V4 are as follows: 



Vi 


= /(") 


/ 


dfL4(n - n A ) 




(28) 




xRe 


l« 


-'t 2 *- 2 *«>(Q(n)+»(n))] 


COs( 27TU ■ 


■n) 


V-2 


= /M 


J 


dQA(h - n.4) 




(29) 




xRe 


[e 


- i(2 *- 2 * (ft)) (Q(n) + it/(n))] 


sin( 27ru • 


ft), 


v 3 


= /M 


J 


d£M(n - n 4 ) 




(30) 




xlm 


[e 


- <(2V, - 2 * Cfl)) (Q(n)+if7(n))] 


COs( 27TU 


•ft), 


V4 


= /M 


J 


dQA(h - n A ) 




(31) 




xlm 




- l(2 ^- 2 * (fl)) (Q(n) + it/(n))] 


sin( 2nu ■ 


ft)- 



With these, we can form two linear combinations Ve and 
Vb such that 

V E = 

ae[ e -«(*»-*)(Vi + i v a )] +iRe[e- t2 ^-^(V 2 + i Vi)], 
V B = 

Imle-^-^iV! +iV 3 )} + ilrn^- 12 ^-^ {V 2 + i Vi)], 
With Eq.[28l[29l|30]and|3Tl we can easily show that 

V E = f(y) J d04(n-fU) 

xRe[e- l2( ' / ' u -* (fl)) (g(n) + i £/(n))]e l2,ru ' fl 
V B = fly) J dnA(n-n A ) 

xIm[e- l2Wu -' I,(fl)) ((3(n) +i(/(n))]e l2,ru,ft . 

We can identify the linear combination Ve with the Vqi 
and Vb with the Vjji , whose polarizer frame is aligned with 
the associated baseline. Hence, the linear combination Ve is 
optimized for E mode measurement and Vb is for B mode 
measurement. Linear combination Ve and Vb enable the 
optimization for a array of multiple feedhorns (N>2). 

7 CONCLUSION 

As shown in this paper, the contribution from E and B mode 
to Vqi and Vjji are variant under the rotation of the polarizer 
frame. In appendix [Bj we have shown analytically that the 
diagonal window functions of Vqi and Vu> are maximized, 
when x axis of the polarizer frame is rotated from the asso- 
ciated baseline with certain angles. The la error on power 
spectra estimation from Vqi and Vu> measurement are also 
variant under the rotation of the polarizer frame. Huge vari- 
ation of la error are shown in the simulated observation 
of the configuration similar to the DASI. The la error on 
B mode power spectrum estimation is minimized, when x 
axis of polarizers is rotated from the baseline by 45° in Vqi 
measurement (when x axis of polarizers is aligned with the 
baseline in Vjji measurement). Simulation shows that mini- 
mum la error on B mode power spectra estimation in Vqi 
or Vjji is 26 per cent of that from Vrl or Vlr, though more 
information are measured in Vrl or Vlr- The simulation 
also shows that the E/B mixing in Vqi or Vyi can be as 



low as 23 per cent of that in Vrl or Vlr- With choice of 
polarizer rotation from the baseline, we can achieve B mode 
sensitivity and E/B separability from Vqi or Vyi measure- 
ment several times better than Vrl or Vlr measurement in 
the same configuration. 

For a array of multiple feedhorns (N>2), there always 
exist certain linear combinations of visibilities, which are 
equivalent to visibilities of the optimal polarizer rotation. 
B mode optimization can be achieved for a feedhorn array 
(N>2) by forming the linear combinations. 

With integration time which makes the noise variance 
equal to the sample varian ce, the estimation error on the 
powe r spectra is minimized l|Bowden et al.l l2004; Par k et ail 
2003). In parallel with the choice for the polarizer rotation 
from the baseline, we can optimize the interferometer for 
B mode with the integration time which makes the sample 
variance of B mode equal to the noise variance. 
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APPENDIX A: CMB POLARIZATION BASIS 
VECTORS AND ANTENNA COORDINATE 

In all-sky analysis, the CMB polarization at the angular 
coordinate (9, <f>) are measured in the local reference frame 
whose axises are (e$, e^, e r ). Let's call this coordinate frame 
'the local CMBP frame' from now on. Consider the polar- 
ization observation of antenna pointing (9a,4>a). A global 
coordinate frame coincides with the antenna coordinate by 
Euler rotations R !/ (Sa) R z (0a). Since a global coordinate 
frame coincides with the local CMBP frame by Euler rota- 
tions R y (6») R z (0), the Euler Rotations R z (j) R M (/3) R z (a) 
coincides the antenna coordinate frame with the local 
CMBP frame, where R z (7) R H (/3) R z (a) R;,(#a) R z (0a) = 
Rj,(0)R z (0). Therefore, the local CMBP frame is in ro- 
tation from the antenna coordinate by the Euler angles 
(a, j3, 7) as follows: 



P 
7 



tan 

COS~~ [C< is I) n is I 

-1 



tan 



a cos(0 — 4>a) — cos 9 sin 9a 
+ sin 9 sin 9a cos(0 — 4>a)] 
sin 9a sin(0 — 0a) 



sin 9 cos 9a + cos 9 sin 9a cos( 



where the Euler angles (a, (3, 7) can be obtained from 
R z ( 7 )R y (/3)R z (a) = R ! ,(6»)R z (0)RJ 1 (0a)R.; 1 (6'a). In 
most CMB polarization experiments, where polarizers are 
attached to the other side of feedhorns, incoming rays go 
through polarizers after feedhorns. After passing through a 
feedhorn, an incoming off-axis ray becomes an on-axis ray. 
Then the local CMBP frame of the ray after the feedhorn 
system is simply in azimuthal rotation a + 7 from the an- 
tenna coordinate. Therefore, $ in Eq. [5] is 



tan 



+ tan 



sin 9 sin(0 — 0a) 



sin 9 cos 9a cos(0 — 0a) — cos 9 sin 9a 
sin 9a sin((/> — 0a) 



— sin 9 cos 9a + cos 9 sin ( 



APPENDIX B: WINDOW FUNCTIONS 

Bl flat sky approximation in small angle limit 

As discussed in Sj4] visibilities with flat sky approximation 
are as follows: 



Vip = 



V RL 



Vlr 



J dVA(u-u') 

x [cos(2(</> - 0»/))B(u') + sin(2(V - U <)) B(u% 
f(u) J dVi(u - u') 

x [- sin(2f> - u ,))E(u') + cos(2(V> - U )) B(u% 
/(«/) /dVi(u-u') 

/(•/) /d 2 u'A(u-u') 

x [e^^-^'^^uO-iSCu'))] • 

For multipole £ > 60, the following relations between the 
flat sky power spectra and the exact power spectra from 
spher ical sky works within one percent error (|White et al.l 
Il999fj : 



(E(u)E*(u')) « Cf B <5(u-u')| ;=2 _, 
<fl(u)B>')> « Cf B «5(u-u')| i=2 _, 
<S(u)B*(u')) = 0. 



(Bl) 
(B2) 
(B3) 



With the correspondence of the power spectra between flat 
sky and spherical sky, it can be easily derived that diago- 
nal elements of E/B window functions and their derivatives 
with respect to the rotation of the polarizer frame, ip, are as 
follows: 

(i) (v Q ,(u)v Q ,(uy), 

W EE (u,u') = 

f{y) u J cos 2 (2i/ ) -20 u /)i(u-u')i(u-u')d0 u /, 

W (U, U ) = 

f{v)u / sin 2 (2-0-20 u ,)i(u-u')i(u-u')d0 u ,. 



f(y) u' J d0 u ,i 2 (u-u') 

x [— 2 sin 4i/> cos 40 u / + 2 cos 4tp sin 40 u /] , 

/>) u' J d^,A 2 (u-u') 

x [2 sin 4-0 cos 40 u / — 2 cos 4?/> sin 40 u /] . 



dip 

dW BB (u,u') 



dip 



(ii) <W'(u)W'(u)*), 
W EE (u,u') = 

/ 3 (z/) u' y sin 2 (2-0 - 20 u ,)i(u - u')A(u - u')d0 u 
W (u, u ) = 



Jaiseung Kim 



f(y)u / cos 2 (2V>-20 u ,)A(u-u')i(u-u')d0 u 



dW EE {u,u') 

dtp 

dW BB {u,u) 

dip 



= f(v)u' I d0 u M 2 (u-"') 

x [2 sin 4ip cos 40 u / — 2 cos 4^/> sin 40 
= />)«' J d0 u , A 2 (u-u') 

x [—2 sin 4?/> cos 40 u / + 2 cos 40 sin ' 



(iii) (V m (u)Vr£(u)*> and {V lr (u)Vlr(u)*), 
W EE (u,u') = W BB (u,u') = 

f{v)u / A(u- u')i(u- u')d0 u ,, 



3W m (u,h') _ dW BB (u,u') 



Of 



From (iii) we can see that the diagonal E and B mode 



window functions in Vrl and Vlr measurement are invari- 



sin 4ip J d0 u /7i 2 (u — u') cos 40 u / 

= cos 4-0 / d0 u / J 4 2 (u — u') sin40 u /]. 



ant under the rotation of the polarizer frame. From (i) and 
(ii)| we can see that the derivatives of the diagonal E/B win- 
dow functions in Vqi and Vv> measurement are zero, when 



(B4) 



ip, which satisfies Eq. IB4I is 

1 _ x / f d0 u , sin40 u ,i 2 (u- u') 

w = — tan 4 ~ 

4 \J d0 u / cos 40 u , A 2 (u - u'] 

A Gaussian beam is a good approximation to many CMB 
experiments and the Fourier transform of a Gaussian beam 
is 

A(u — u') = exp[— - 



(B5) 



2 

2 , _./2 



= exp[- 



U + U — 2UU COs(0 u / — 0u)]<7 



where a = 0.4245 FWHM. With a Gaussian beam, the ar- 
gument of tan" 1 in Eg. IB5I is 

J 2 ^ d0 u / sin40 u / exp[2ww' cos(0 u / — 0u)<?" 2 ] 
J 2n d0 u / cos40 u / exp[2mt' cos(0 u / — u )<t 2 ] 



d0 sin 4(0 + U ) exp[2uu' cos(cp)a 2 



2w — t 
2tt 



d0 cos 4(0 + U ) exp[2im' cos(0)<r 2 



J * dtp sin 4(0 + U ) exp[2««' cos(0)er 2 ] 
J 2 ^ d(p cos 4(0 + U ) exp[2tt-u' cos(0)<r 2 ] 
sin40 u J" 2lr d0 cos 40exp[2mt' cos(0)<r 2 ] 
cos 40 u J 2n d0 cos40exp[2mi' cos(0)<j 2 ] 
tan40 u 



From the third line to the fourth line in Eq. IB6I 
J n d0sin 40 exp[2uu' cos(0)<j 2 ] = was used. By plugging 



(B6) 



Eq. [B6]into Eq. [B5] we get 

ip = — tan -1 (tan40 u ) 

727T 

= 0u + -r 



(n 



where U is the orientation of the baseline. The diagonal 
element of E and B mode window functions of Vqi and Vw 
is maximized or minimized, when x axis of the polarizer is in 
rotation from the baseline by -90°, -45°, 0°, 45°, 90°. Since 
the second derivative of diagonal element of E mode window 
function has the opposite sign of that of B mode window 
function, the diagonal element of B mode window function 
is minimized at the polarizer rotation which maximizes the 
diagonal element of E mode window function, and vice versa. 



B2 spherical sky 

Visibilities from spherical sky are as follows: 



V Q , = / dtM(n,fu)e'< 2 ™- fl ) 

,,r — i(2i/> — 2*(n)) / . • \ v 

x[e (a E ,i m + ia B ,im) %Yfr, 



+e 



i(2-0-2#(fl)) 



(a E ,l m ?■ @>B ,lm) — 2 i/mj 



Vu> = j/W / d^(n,n,)e ,(2 "- 9) 



(B7) 



(B8) 



-i(2t/>-2$(A)) 



Pin)) / | \ -tr 

(0,E,lm + lClB,lm) 211m 

-e y v (a E ,i m - ia B ,im) -nYi m \ 

Vrl = -f(y) J d04(n,n A ) (B9) 

w I i • \ v i(27ru-ft-2j/i+2$(ft)) 

X {dE.lm + I ClB,l m ) 2 J- im£ ; 

Vlr = -f(u) J d04(fi,A A ) (BIO) 

I ■ \ V i(27ru-ft+2;/j-2*(n)) 



With Eq. |B7] EU [B9l and iBlOl visibilities from spherical 
sky can be expressed as follows: 

V Q r(n,u) = y^(e~' 2V, _R; m + e 2i) Li m )a E ,im 

l,m 

+i {e~ l2,p Ri m ~ e 2 ^ L im ) a B ,i m , (Bll) 
V(7'(n,u) = ^ -~i(e~ l2i ' R lm - e n 'Li m )a E ,im 

l,m 

+ (e-' 2 ^R lm + e l2 ^L lm )a B:tm , (B12) 
V RL (n,u) = 2J 2 ( a E,im + ia B ,im)e~ l2 ' 1 ' Rim, (B13) 

i,m 

Vtn(n,u) = 2J 2 (a B j m + ia s , !m )e l2,/ 'L; m , (B14) 

Z ,m 

where 



-Kzm — 



2 

/O) 



dQ.A{n, TLA) 2Y tm e 1 1 v 
dUA(n, n A ) -a^ime k 1 k 



and fiA indicates the direction of antenna pointing. 
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Spin ±2 sphe rical harmonics have the following form 
i|Zaldarriagalll99SF) . 



2 F im (n = 



21 + 1 

47T 

21 + 1 

47T 



(F ldm (e) + F 2 , lm (0)) 



{Fi,i m {6)-F % i m {e))e 



inxp 



where Fijm and F 2 ,i m can be computed in terms of Legendre 
functions as follows ( Kamion kowski et al.|[l997h : 



F Um (9) 



/(/ - 2)!(7 - m)' , , cosfl „ m , 

-( l —^r + h(i-i))pr(cose)], 



F2,lm{6) — 2 



(I - 2)1(1 - m)\ m 
(/ + 2)!(Z + m)!sin T i 

(I- 1) cos 6Pr (cos 6)}. 



■[(l + m)P™x(cost 



y^Me* 4-0 L lm R* m - He m ' R lm L*„ 



(iii) (V RL (u)V RL (u)*), 

W EE = W BB = 4R lm Ri 



dW t EE 
dip 



dijj 



= o. 



(iv) {v LR (u)v LR (uy}, 

W EE = W BB = AL lm L^ m . 



dw EE 

dip 



dWp 



dip 



= 0. 



The covariance properties of the E and B mode are given by 

( a E,lmCL E ,l l m') = C EE 5 n i5 mrn > , 
(o.B,lm,a Bt li m i) 
(o,E,lmdB,l'm') = 0. 



Cf B &ll'8 mrn i 



With these covariance properties, Eq. IBTTI IB121 [B13l and 
IB141 it can be easily shown that diagonal elements of E/B 
window functions and their derivatives with respect to the 
rotation of the polarizer frame, ip, are as follows: 

(i) {VcyfufiVtyW), 

Wj = ^ ' RimRj,im + Li m Lj tm 



Wf 



dW L 



^ ^ RlmRlm + Ll m Li m 



-Jim -rt-Z r, 



-i4ip 



dip 

dw EB 



22 *4e l4,/ 'LlmR'lr, 



i4e ' Ri m L* m , 



dip 

m 

(ii) (V v ,(u)Vu,(ur) 



-i^Li m Rj m + He~^R lm L\„ 



Wt 



^ ' RlmRlm + Ll 

m 

— e ^Im-tilm ~~ e K'lm-L'lmi 
^ RlmRlm + Ll m L t „ 



+ e Lilmttlm + e ttlmlJlr, 



From [(iilJI and [(iv)] we can see that the diagonal E and 
B mode window functions in Vbl and Vlr measurement are 



invariant under the rotation of the polarizer frame. From (i) 



and |(ii)| we can see that the derivatives of the diagonal E 
and B mode window functions in Vqi and Vu> measurement 
are zero, when 



f r»* —Hip ry j * 

6 ^i^m-ttijm ^ ^i,lm^i Im 



0. (B15) 

Since the left side of Eq. IB15I is 

-i2Im[e"' 4,l, fl, i i m L* lm ] 

= -i2(- sin&ipRe[Ri t i m L it i m *] + cos4i/dm [Ri,i m L iilm *]) , 

the following ip satisfies Eq. IB15I 



1 _i Im [RijmLij 
ip = — tan — — r-^- 2 — ^ — - 

4 Re [RiAmLi Am 



(B16) 



The diagonal elements of window function are in- 
dependent of the choice of the reference coordinate 
ijWhite and Srednickil Il995h . So we can choose antenna 
pointing as z axis without loss of generality. In the refer- 
ence frame of our choice, < l > (n) is reduced to azimuthal an- 
gle 4>. In most of CMB interferometer experiments, primary 
beampattern are azimuthally symmetric and baselines are 
coplanar. Then Ri m and L; m are 



R, 



fiy) 



x 2 Yi 



2Yim{V,(p)e 



d(0) s\n6 A(6) 
<i Jo 

i(2iru sin 9 cos(</>— U ) — 2i/)+20) 



f(v) 1 21 + 1 
x / d(6) sine A(6)[F 1 ,im(0) + F 2 jm(6)} 

>0 

f2n 



X 



o i(2iru sin B cos(<^ — </> u ) — 2^ + (m + 2)<^>) 



if (m + 2) i(m + 2)</> u r\ 

c 2v- , Zm, 



(B17) 



9V 



J i,lm 



fiy) 



A(6) sinO A(6) I d<j> 
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X- 2 Y lm (e, <j>)e i( - 2lru sin 9 c^-M+^+W) ) 
f(y) 121 + 1 







4tt 

■ / d(9)sm8A(9)[F Um (0) - F Vm (0)] 

. J ^,i(2TTUsili6co$(<t>-4> l i)+2'ij>+(m—2)<t>) 


e^ (m - 2) e l(m - 2)</> "- 2 e im , (B18) 



where 



.(-),,„ = -f{v)\l ( 2l + 1 > [ d9sm(e)A(e) 



(J 



x J m+2 (27rusin(0))(F 1 , im (6>) + F 2 , lm (9)), 



,(->,,„ = -/(i/) W ( a + 1 > ,r / d0sin(0)A(0) 



x J m _ 2 (27rasin(0))(Fi, im (0) - F 2>im (0)). 



By plugging Eg . IB171 and [BUI into Eg. [B16l we get 

1 . _i Im [Ri,imLii m I 
ip — — tan 



I Re \Ri imLi i 

I rn 



i t _x im[e i4 ^] 2 e im _ 2 e ; . 



— tan 



4 Re[e <4 *u] 2 e Jm _ 2 e im 

= 0u + ^. (n = ... ,-2, -1,0,1,2,...) 

^ = 0u + ^f- maximizes or minimizes the diagonal element 
of E and B mode window functions of Vqi and Vw , which 
is consistent with the result obtained with flat sky approxi- 
mation. 



